AAPM WGDCAB Report 372: A joint AAPM, ESTRO, ABG, and ABS report on commissioning of model-based dose calculation algorithms in brachytherapy.

The introduction of model-based dose calculation algorithms (MBDCAs) in brachytherapy provides an opportunity for a more accurate dose calculation and opens the possibility for novel, innovative treatment modalities. The joint AAPM, ESTRO, and ABG Task Group 186 (TG-186) report provided guidance to early adopters. However, the commissioning aspect of these algorithms was described only in general terms with no quantitative goals. This report, from the Working Group on Model-Based Dose Calculation Algorithms in Brachytherapy, introduced a field-tested approach to MBDCA commissioning. It is based on a set of well-characterized test cases for which reference Monte Carlo (MC) and vendor-specific MBDCA dose distributions are available in a Digital Imaging and Communications in Medicine-Radiotherapy (DICOM-RT) format to the clinical users. The key elements of the TG-186 commissioning workflow are now described in detail, and quantitative goals are provided. This approach leverages the well-known Brachytherapy Source Registry jointly managed by the AAPM and the Imaging and Radiation Oncology Core (IROC) Houston Quality Assurance Center (with associated links at ESTRO) to provide open access to test cases as well as step-by-step user guides. While the current report is limited to the two most widely commercially available MBDCAs and only for 192 Ir-based afterloading brachytherapy at this time, this report establishes a general framework that can easily be extended to other brachytherapy MBDCAs and brachytherapy sources. The AAPM, ESTRO, ABG, and ABS recommend that clinical medical physicists implement the workflow presented in this report to validate both the basic and the advanced dose calculation features of their commercial MBDCAs. Recommendations are also given to vendors to integrate advanced analysis tools into their brachytherapy treatment planning system to facilitate extensive dose comparisons. The use of the test cases for research and educational purposes is further encouraged.

A MC-based anthropomorphic test case for commissioning model-based dose calculation in interstitial breast 192-Ir HDR brachytherapy.

PURPOSE: To provide the first clinical test case for commissioning of 192 Ir brachytherapy model-based dose calculation algorithms (MBDCAs) according to the AAPM TG-186 report workflow. ACQUISITION AND VALIDATION METHODS: A computational patient phantom model was generated from a clinical multi-catheter 192 Ir HDR breast brachytherapy case. Regions of interest (ROIs) were contoured and digitized on the patient CT images and the model was written to a series of DICOM CT images using MATLAB. The model was imported into two commercial treatment planning systems (TPSs) currently incorporating an MBDCA. Identical treatment plans were prepared using a generic 192 Ir HDR source and the TG-43-based algorithm of each TPS. This was followed by dose to medium in medium calculations using the MBDCA option of each TPS. Monte Carlo (MC) simulation was performed in the model using three different codes and information parsed from the treatment plan exported in DICOM radiation therapy (RT) format. Results were found to agree within statistical uncertainty and the dataset with the lowest uncertainty was assigned as the reference MC dose distribution. DATA FORMAT AND USAGE NOTES: The dataset is available online at http://irochouston.mdanderson.org/rpc/BrachySeeds/BrachySeeds/index.html,https://doi.org/10.52519/00005. Files include the treatment plan for each TPS in DICOM RT format, reference MC dose data in RT Dose format, as well as a guide for database users and all files necessary to repeat the MC simulations. POTENTIAL APPLICATIONS: The dataset facilitates the commissioning of brachytherapy MBDCAs using TPS embedded tools and establishes a methodology for the development of future clinical test cases. It is also useful to non-MBDCA adopters for intercomparing MBDCAs and exploring their benefits and limitations, as well as to brachytherapy researchers in need of a dosimetric and/or a DICOM RT information parsing benchmark. Limitations include specificity in terms of radionuclide, source model, clinical scenario, and MBDCA version used for its preparation.

Dosimetric accuracy of the Convolution algorithm for Leksell Gamma Plan radiosurgery treatment planning: Evaluation in the presence of clinically relevant inhomogeneities.

PURPOSE: The Leksell Gamma Plan Convolution algorithm (LGP-Convolution) has not been widely adopted. This mainly stems from the higher calculated beam-on times relative to the standard ray tracing-based LGP-TMR10 dose calculation algorithm. This study aims to evaluate the accuracy of the LGP-Convolution in scenarios where the treated lesions are in the vicinity of or encompassed by bone and/or air inhomogeneities. METHODS: The solid water dosimetry phantom provided by the vendor was modified to include bone and air inhomogeneities. Two treatment planning scenarios were investigated involving a single shot and multiple shots, respectively. Treatment planning and dose prescription were performed using the LGP-Convolution algorithm. Triple channel film dosimetry was performed using GafChromic EBT3 films calibrated in terms of absorbed dose to water in a 60 Co beam. Monte Carlo (MC) simulation dosimetry was also performed in the inhomogeneous experimental geometry using the EGSnrc MC platform and a previously validated sector-based phase-space source model. MC simulations were also employed to determine correction factors required for converting EBT3 measurements at points within the bone and air inhomogeneities from dose-to-water values to the corresponding dose to medium values. RESULTS AND CONCLUSIONS: EBT3 dose to medium correction factors ranged with field size (4, 8, or 16 mm) within 0.941-0.946 for bone and 0.745-0.749 for air inhomogeneities. An excellent agreement was found between the LGP-Convolution calculations with corresponding EBT3 and MC dose to medium results at all measurement points, except those located inside the air inhomogeneity. The latter is of no clinical importance and excluding them yielded gamma index passing rates of nearly 100% for 3% local dose difference and 1 mm distance-to-agreement criteria. The excellent agreement observed between LGP-Convolution calculations and film as well as MC results of dose to medium indicates that the latter is the quantity reported by the LGP-Convolution.

Dosimetry in 1.5 T MR-Linacs: Monte Carlo determination of magnetic field correction factors and investigation of the air gap effect.

BACKGROUND: In Magnetic Resonance-Linac (MR-Linac) dosimetry formalisms, a new correction factor, kB,Q , has been introduced to account for corresponding changes to detector readings under the beam quality, Q, and the presence of magnetic field, B. PURPOSE: This study aims to develop and implement a Monte Carlo (MC)-based framework for the determination of kB,Q correction factors for a series of ionization chambers utilized for dosimetry protocols and dosimetric quality assurance checks in clinical 1.5 T MR-Linacs. Their dependencies on irradiation setup conditions are also investigated. Moreover, to evaluate the suitability of solid phantoms for dosimetry checks and end-to-end tests, changes to the detector readings due to the presence of small asymmetrical air gaps around the detector's tip are quantified. METHODS: Phase space files for three irradiation fields of the ELEKTA Unity 1.5 T/7 MV flattening-filter-free MR-Linac were provided by the manufacturer and used as source models throughout this study. Twelve ionization chambers (three farmer-type and nine small-cavity detectors, from three manufacturers) were modeled (including their dead volume) using the EGSnrc MC code package. kB,Q values were calculated for the 10 × 10 cm2 irradiation field and for four cardinal orientations of the detectors' axes with respect to the 1.5 T magnetic field. Potential dependencies of kB,Q values with respect to field size, depth, and phantom material were investigated by performing additional simulations. Changes to the detectors' readings due to the presence of small asymmetrical air gaps (0.1 up to 1 mm) around the chambers' sensitive volume in an RW3 solid phantom were quantified for three small-cavity chambers and two orientations. RESULTS: For both parallel (to the magnetic field) orientations, kB,Q values were found close to unity. The maximum correction needed was 1.1%. For each detector studied, the kB,Q values calculated for the two parallel orientations agreed within uncertainties. Larger corrections (up to 5%) were calculated when the detectors were oriented perpendicularly to the magnetic field. Results were compared with corresponding ones found in the literature, wherever available. No considerable dependence of kB,Q with respect to field size (down to 3 × 3 cm2 ), depth, or phantom material was noticed, for the detectors investigated. As compared to the perpendicular one, in the parallel to the magnetic field orientation, the air gap effect is minimized but is still considerable even for the smallest air gap considered (0.1 mm). CONCLUSION: For the 10 × 10 cm2 field, magnetic field correction factors for 12 ionization chambers and four orientations were determined. For each detector, the kB,Q value may be also applied for dosimetry procedures under different irradiation parameters provided that the orientation is taken into account. Moreover, if solid phantoms are used, even the smallest asymmetrical air gap may still bias small-cavity chamber response. This work substantially expands the availability and applicability of kB,Q correction factors that are detector- and orientation-specific, enabling more options in MR-Linac dosimetry checks, end-to-end tests, and quality assurance protocols.

Spatial distribution of the imaging dose and characterization of the scatter radiation contribution in CyberKnife radiosurgery.

PURPOSE: The imaging dose for intra- and extra-cranial CyberKnife radiosurgery applications was calculated and the scattered radiation reaching the digital detectors was quantified and analyzed with regard to its origin. METHODS: The image guidance subsystem of the CyberKnife was modeled based on vendor-provided information. The emitted X-ray energy spectrum for 120 kV was estimated using the SpekPy software tool. Monte Carlo (MC) image acquisition simulations were performed to calculate the total, primary and scattered photon fluences reaching each detector as a function of the imaged object dimensions. MC calculations of the imaging dose were performed for intra- and extra-cranial applications assuming 120 kV and 10 mAs acquisition settings. RESULTS: The amount of scattered radiation reaching each detector was found to depend on the dimensions of the imaged anatomical region, contributing more than 40 % to the total photon fluence for regions more than 20 cm thick. More than 20 % of this scattered radiation originates from the contralateral imaging field. A maximum organ dose of 1.5 mGy at the nasal bones and an average dose of 0.37 mGy to the eye lenses per image pair acquisition was calculated for head applications. An entrance imaging dose of 0.4 mGy was calculated for extracranial applications. CONCLUSIONS: Scattered radiation reaching each detector in the skull and spine tracking applications can be reduced by acquiring the pair of radiographs sequentially instead of simultaneously. A dose of 3.7 cGy to the eye lenses is estimated assuming 100 image pair exposures required for treatment completion.

On the potential of 2D ion chamber arrays for high-dose rate remote afterloading brachytherapy quality assurance.

Objective. To investigate the potential of 2D ion chamber arrays to serve as a standalone tool for the verification of source strength, positioning and dwell time, within the framework of192Ir high-dose rate brachytherapy device quality assurance (QA).Approach.A commercially available ion chamber array was used. Fitting of a 2D Lorentzian peak function to experimental data from a multiple source dwell position irradiation on a frame-by-frame basis, facilitated tracking of the source center orthogonal projection on the array plane. For source air kerma strength verification, Monte Carlo simulation was employed to obtain a chamber array- and source-specific correction factor of calibration with a 6 MV photon beam. This factor converted the signal measured by each ion chamber element to air kerma in free space. A source positioning correction was also applied to lift potential geometry mismatch between experiment and Monte Carlo simulation.Main results.Spatial and temporal accuracy of source movement was verified within 0.5 mm and 0.02 s, respectively, in compliance with the test endpoints recommended by international professional societies. The source air kerma strength was verified experimentally within method uncertainties estimated as 1.44% (k = 1). The source positioning correction method employed did not introduce bias to experimental results of irradiations where source positioning was accurate. Development of a custom jig attachable to the chamber array for accurate and reproducible experimental set up would improve testing accuracy and obviate the need for source positioning correction in air kerma strength verification.Significance.Delivery of a single irradiation plan, optimized based on results of this work, to a 2D ion chamber array can be used for concurrent testing of source position, dwell time and air kerma strength, and the procedure can be expedited through automation. Chamber arrays merit further study in treatment planning QA and real time,in vivodose verification.

A comparative assessment of inhomogeneity and finite patient dimension effects in 60Co and 192Ir high-dose-rate brachytherapy.

PURPOSE: To perform a comparative study of heterogeneities and finite patient dimension effects in 60Co and 192Ir high-dose-rate (HDR) brachytherapy. MATERIAL AND METHODS: Clinically equivalent plans were prepared for 19 cases (8 breast, 5 esophagus, 6 gynecologic) using the Ir2.A85-2 and the Co0.A86 HDR sources, with a TG-43 based treatment planning system (TPS). Phase space files were obtained for the two source designs using MCNP6, and validated through comparison to a single source dosimetry results in the literature. Dose to water, taking into account the patient specific anatomy and materials (Dw,m), was calculated for all plans using MCNP6, with input files prepared using the BrachyGuide software tool to analyze information from DICOM RT plan exports. RESULTS: A general TG-43 dose overestimation was observed, except for the lungs, with a greater magnitude for 192Ir. The distribution of percentage differences between TG-43 and Monte Carlo (MC) in dose volume histogram (DVH) indices for the planning target volume (PTV) presented small median values (about 2%) for both 60Co and 192Ir, with a greater dispersion for 192Ir. Regarding the organs at risk (OARs), median percentage differences for breast V50% were 3% (5%) for 60Co (192Ir). Differences in median skin D2cc were found comparable, with a larger dispersion for 192Ir, and the same applied to the lung D10cc and the aorta D2cc. TG-43 overestimates D2cc for the rectum and the sigmoid, with median differences from MC within 2% and a greater dispersion for 192Ir. For the bladder, the median of the difference is greater for 60Co (~2%) than for 192Ir (~0.75%), demonstrating however a greater dispersion again for 192Ir. CONCLUSIONS: The magnitude of differences observed between TG-43 based and MC dosimetry and their smaller dispersion relative to 192Ir, suggest that 60Co HDR sources are more amenable to the TG-43 assumptions in clinical treatment planning dosimetry.

On the use of a novel Ferrous Xylenol-orange gelatin dosimeter for HDR brachytherapy commissioning and quality assurance testing.

PURPOSE: To evaluate a commercially available Ferrous-Xylenol Orange-Gel (FXG) dosimeter (TrueView™) coupled with Optical-Computed Tomography (OCT) read out, for 3D dose verification in an Ir-192 superficial brachytherapy application. METHODS: Two identical polyethylene containers filled with gel from the same batch were used. One was irradiated with an 18 MeV electron field to examine the dose-response linearity and obtain a calibration curve. A flap surface applicator was attached to the other to simulate treatment of a skin lesion. The dose distribution in the experimental set up was calculated with the TG-43 and the model based dose calculation (MBCA) algorithms of a commercial treatment planning system (TPS), as well as Monte Carlo (MC) simulation using the MCNP code. Measured and calculated dose distributions were spatially registered and compared. RESULTS: Apart from a region close to the container's neck, where gel measurements exhibited an over-response relative to MC calculations (probably due to stray light perturbation), an excellent agreement was observed between measurements and simulations. More than 97% of points within the 10% isodose line (80 cGy) met the gamma index criteria established from uncertainty analysis (5%/2 mm). The corresponding passing rates for the comparison of experiment to calculations using the TG-43 and MBDCA options of the TPS were 57% and 92%, respectively. CONCLUSION: TrueView™ is suitable for the quality assurance of demanding radiotherapy applications. Experimental results of this work confirm the advantage of the studied MBDCA over TG-43, expected from the improved account of scatter radiation in the treatment geometry.

A generic TG-186 shielded applicator for commissioning model-based dose calculation algorithms for high-dose-rate 192 Ir brachytherapy.

PURPOSE: A joint working group was created by the American Association of Physicists in Medicine (AAPM), the European Society for Radiotherapy and Oncology (ESTRO), and the Australasian Brachytherapy Group (ABG) with the charge, among others, to develop a set of well-defined test case plans and perform calculations and comparisons with model-based dose calculation algorithms (MBDCAs). Its main goal is to facilitate a smooth transition from the AAPM Task Group No. 43 (TG-43) dose calculation formalism, widely being used in clinical practice for brachytherapy, to the one proposed by Task Group No. 186 (TG-186) for MBDCAs. To do so, in this work a hypothetical, generic high-dose rate (HDR) 192 Ir shielded applicator has been designed and benchmarked. METHODS: A generic HDR 192 Ir shielded applicator was designed based on three commercially available gynecological applicators as well as a virtual cubic water phantom that can be imported into any DICOM-RT compatible treatment planning system (TPS). The absorbed dose distribution around the applicator with the TG-186 192 Ir source located at one dwell position at its center was computed using two commercial TPSs incorporating MBDCAs (Oncentra® Brachy with Advanced Collapsed-cone Engine, ACE™, and BrachyVision ACUROS™) and state-of-the-art Monte Carlo (MC) codes, including ALGEBRA, BrachyDose, egs_brachy, Geant4, MCNP6, and Penelope2008. TPS-based volumetric dose distributions for the previously reported "source centered in water" and "source displaced" test cases, and the new "source centered in applicator" test case, were analyzed here using the MCNP6 dose distribution as a reference. Volumetric dose comparisons of TPS results against results for the other MC codes were also performed. Distributions of local and global dose difference ratios are reported. RESULTS: The local dose differences among MC codes are comparable to the statistical uncertainties of the reference datasets for the "source centered in water" and "source displaced" test cases and for the clinically relevant part of the unshielded volume in the "source centered in applicator" case. Larger local differences appear in the shielded volume or at large distances. Considering clinically relevant regions, global dose differences are smaller than the local ones. The most disadvantageous case for the MBDCAs is the one including the shielded applicator. In this case, ACUROS agrees with MC within [-4.2%, +4.2%] for the majority of voxels (95%) while presenting dose differences within [-0.12%, +0.12%] of the dose at a clinically relevant reference point. For ACE, 95% of the total volume presents differences with respect to MC in the range [-1.7%, +0.4%] of the dose at the reference point. CONCLUSIONS: The combination of the generic source and generic shielded applicator, together with the previously developed test cases and reference datasets (available in the Brachytherapy Source Registry), lay a solid foundation in supporting uniform commissioning procedures and direct comparisons among treatment planning systems for HDR 192 Ir brachytherapy.

Supplement 2 for the 2004 update of the AAPM Task Group No. 43 Report: Joint recommendations by the AAPM and GEC-ESTRO.

Since the publication of the 2004 update to the American Association of Physicists in Medicine (AAPM) Task Group No. 43 Report (TG-43U1) and its 2007 supplement (TG-43U1S1), several new low-energy photon-emitting brachytherapy sources have become available. Many of these sources have satisfied the AAPM prerequisites for routine clinical purposes and are posted on the Brachytherapy Source Registry managed jointly by the AAPM and the Imaging and Radiation Oncology Core Houston Quality Assurance Center (IROC Houston). Given increasingly closer interactions among physicists in North America and Europe, the AAPM and the Groupe Européen de Curiethérapie-European Society for Radiotherapy & Oncology (GEC-ESTRO) have prepared another supplement containing recommended brachytherapy dosimetry parameters for eleven low-energy photon-emitting brachytherapy sources. The current report presents consensus datasets approved by the AAPM and GEC-ESTRO. The following sources are included: 125 I sources (BEBIG model I25.S17, BEBIG model I25.S17plus, BEBIG model I25.S18, Elekta model 130.002, Oncura model 9011, and Theragenics model AgX100); 103 Pd sources (CivaTech Oncology model CS10, IBt model 1031L, IBt model 1032P, and IsoAid model IAPd-103A); and 131 Cs (IsoRay Medical model CS-1 Rev2). Observations are included on the behavior of these dosimetry parameters as a function of radionuclide. Recommendations are presented on the selection of dosimetry parameters, such as from societal reports issuing consensus datasets (e.g., TG-43U1, AAPM Report #229), the joint AAPM/IROC Houston Registry, the GEC-ESTRO website, the Carleton University website, and those included in software releases from vendors of treatment planning systems. Aspects such as timeliness, maintenance, and rigor of these resources are discussed. Links to reference data are provided for radionuclides (radiation spectra and half-lives) and dose scoring materials (compositions and mass densities). The recent literature is examined on photon energy response corrections for thermoluminescent dosimetry of low-energy photon-emitting brachytherapy sources. Depending upon the dosimetry parameters currently used by individual physicists, use of these recommended consensus datasets may result in changes to patient dose calculations. These changes must be carefully evaluated and reviewed with the radiation oncologist prior to their implementation.

Time resolved dose rate distributions in brachytherapy.

PURPOSE: To investigate the biological significance of introducing time-resolved dose rate distributions (TR-DRD) in brachytherapy. MATERIALS AND METHODS: The treatment plan of a head and neck patient treated with pulsed-dose-rate (PDR) brachytherapy was considered. The TR-DRD was calculated on the basis of a Monte Carlo generated single source dose rate matrix taking into account the dose rate per source dwell position. Biologically Effective Dose (BED) was obtained considering either the mean dose rate per pulse (analytical method) or the TR-DRD (numerical method). Corresponding Tumor Control Probabilities (TCP) were calculated and compared for various PDR schemes and repair half-times from the literature. The dose of the biologically equivalent high-dose-rate (HDR) treatment schedule was also evaluated. RESULTS: The analytical method presents an overall BED underestimation (up to 2%) relative to TR-DRD results. This is associated with an analytical-based TCP underestimation which increases with dose/pulse, pulse duration and period time and decreases with total dose. The half-time of repair seems to have the largest impact on the TCP calculations, with significant differences (up to 39.1%) corresponding to the shorter repair half-times. Regarding the equivalent HDR treatment schedule, the analytical method resulted to a HDR isoeffective dose underestimation lower than 2.2% and thus does not warrant any change in the derivation of the equivalent HDR scheme. CONCLUSION: TR-DRD data should be taken into account for PDR biological effectiveness estimations, especially for short tissue repair half-times. This does not appear however to influence dose prescription of the equivalent HDR treatment schedule for mobile tongue carcinoma.

On the experimental validation of model-based dose calculation algorithms for 192Ir HDR brachytherapy treatment planning.

There is an acknowledged need for the design and implementation of physical phantoms appropriate for the experimental validation of model-based dose calculation algorithms (MBDCA) introduced recently in 192Ir brachytherapy treatment planning systems (TPS), and this work investigates whether it can be met. A PMMA phantom was prepared to accommodate material inhomogeneities (air and Teflon), four plastic brachytherapy catheters, as well as 84 LiF TLD dosimeters (MTS-100M 1 × 1 × 1 mm3 microcubes), two radiochromic films (Gafchromic EBT3) and a plastic 3D dosimeter (PRESAGE). An irradiation plan consisting of 53 source dwell positions was prepared on phantom CT images using a commercially available TPS and taking into account the calibration dose range of each detector. Irradiation was performed using an 192Ir high dose rate (HDR) source. Dose to medium in medium, [Formula: see text], was calculated using the MBDCA option of the same TPS as well as Monte Carlo (MC) simulation with the MCNP code and a benchmarked methodology. Measured and calculated dose distributions were spatially registered and compared. The total standard (k = 1) spatial uncertainties for TLD, film and PRESAGE were: 0.71, 1.58 and 2.55 mm. Corresponding percentage total dosimetric uncertainties were: 5.4-6.4, 2.5-6.4 and 4.85, owing mainly to the absorbed dose sensitivity correction and the relative energy dependence correction (position dependent) for TLD, the film sensitivity calibration (dose dependent) and the dependencies of PRESAGE sensitivity. Results imply a LiF over-response due to a relative intrinsic energy dependence between 192Ir and megavoltage calibration energies, and a dose rate dependence of PRESAGE sensitivity at low dose rates (<1 Gy min-1). Calculations were experimentally validated within uncertainties except for MBDCA results for points in the phantom periphery and dose levels <20%. Experimental MBDCA validation is laborious, yet feasible. Further work is required for the full characterization of dosimeter response for 192Ir and the reduction of experimental uncertainties.

On the impact of improved dosimetric accuracy on head and neck high dose rate brachytherapy.

PURPOSE: To study the effect of finite patient dimensions and tissue heterogeneities in head and neck high dose rate brachytherapy. METHODS AND MATERIALS: The current practice of TG-43 dosimetry was compared to patient specific dosimetry obtained using Monte Carlo simulation for a sample of 22 patient plans. The dose distributions were compared in terms of percentage dose differences as well as differences in dose volume histogram and radiobiological indices for the target and organs at risk (mandible, parotids, skin, and spinal cord). RESULTS: Noticeable percentage differences exist between TG-43 and patient specific dosimetry, mainly at low dose points. Expressed as fractions of the planning aim dose, percentage differences are within 2% with a general TG-43 overestimation except for the spine. These differences are consistent resulting in statistically significant differences of dose volume histogram and radiobiology indices. Absolute differences of these indices are however small to warrant clinical importance in terms of tumor control or complication probabilities. CONCLUSIONS: The introduction of dosimetry methods characterized by improved accuracy is a valuable advancement. It does not appear however to influence dose prescription or call for amendment of clinical recommendations for the mobile tongue, base of tongue, and floor of mouth patient cohort of this study.

A user-oriented procedure for the commissioning and quality assurance testing of treatment planning system dosimetry in high-dose-rate brachytherapy.

PURPOSE: To develop a user-oriented procedure for testing treatment planning system (TPS) dosimetry in high-dose-rate brachytherapy, with particular focus to TPSs using model-based dose calculation algorithms (MBDCAs). METHODS AND MATERIALS: Identical plans were prepared for three computational models using two commercially available systems and the same (192)Ir source. Reference dose distributions were obtained for each plan using the MCNP v.6.1 Monte Carlo (MC) simulation code with input files prepared via automatic parsing of plan information using a custom software tool. The same tool was used for the comparison of reference dose distributions with corresponding MBDCA exports. RESULTS: The single source test case yielded differences due to the MBDCA spatial discretization settings. These affect points at relatively increased distance from the source, and they are abated in test cases with multiple source dwells. Differences beyond MC Type A uncertainty were also observed very close to the source(s), close to the test geometry boundaries, and within heterogeneities. Both MBDCAs studied were found equivalent to MC within 5 cm from the target volume for a clinical breast brachytherapy test case. These are in agreement with previous findings of MBDCA benchmarking in the literature. CONCLUSIONS: The data and the tools presented in this work, that are freely available via the web, can serve as a benchmark for advanced clinical users developing their own tests, a complete commissioning procedure for new adopters of currently available TPSs using MBDCAs, a quality assurance testing tool for future updates of already installed TPSs, or as an admission prerequisite in multicentric clinical trials.

A generic high-dose rate (192)Ir brachytherapy source for evaluation of model-based dose calculations beyond the TG-43 formalism.

PURPOSE: In order to facilitate a smooth transition for brachytherapy dose calculations from the American Association of Physicists in Medicine (AAPM) Task Group No. 43 (TG-43) formalism to model-based dose calculation algorithms (MBDCAs), treatment planning systems (TPSs) using a MBDCA require a set of well-defined test case plans characterized by Monte Carlo (MC) methods. This also permits direct dose comparison to TG-43 reference data. Such test case plans should be made available for use in the software commissioning process performed by clinical end users. To this end, a hypothetical, generic high-dose rate (HDR) (192)Ir source and a virtual water phantom were designed, which can be imported into a TPS. METHODS: A hypothetical, generic HDR (192)Ir source was designed based on commercially available sources as well as a virtual, cubic water phantom that can be imported into any TPS in DICOM format. The dose distribution of the generic (192)Ir source when placed at the center of the cubic phantom, and away from the center under altered scatter conditions, was evaluated using two commercial MBDCAs [Oncentra(®) Brachy with advanced collapsed-cone engine (ACE) and BrachyVision ACUROS™ ]. Dose comparisons were performed using state-of-the-art MC codes for radiation transport, including ALGEBRA, BrachyDose, GEANT4, MCNP5, MCNP6, and PENELOPE2008. The methodologies adhered to recommendations in the AAPM TG-229 report on high-energy brachytherapy source dosimetry. TG-43 dosimetry parameters, an along-away dose-rate table, and primary and scatter separated (PSS) data were obtained. The virtual water phantom of (201)(3) voxels (1 mm sides) was used to evaluate the calculated dose distributions. Two test case plans involving a single position of the generic HDR (192)Ir source in this phantom were prepared: (i) source centered in the phantom and (ii) source displaced 7 cm laterally from the center. Datasets were independently produced by different investigators. MC results were then compared against dose calculated using TG-43 and MBDCA methods. RESULTS: TG-43 and PSS datasets were generated for the generic source, the PSS data for use with the ace algorithm. The dose-rate constant values obtained from seven MC simulations, performed independently using different codes, were in excellent agreement, yielding an average of 1.1109 ± 0.0004 cGy/(h U) (k = 1, Type A uncertainty). MC calculated dose-rate distributions for the two plans were also found to be in excellent agreement, with differences within type A uncertainties. Differences between commercial MBDCA and MC results were test, position, and calculation parameter dependent. On average, however, these differences were within 1% for ACUROS and 2% for ace at clinically relevant distances. CONCLUSIONS: A hypothetical, generic HDR (192)Ir source was designed and implemented in two commercially available TPSs employing different MBDCAs. Reference dose distributions for this source were benchmarked and used for the evaluation of MBDCA calculations employing a virtual, cubic water phantom in the form of a CT DICOM image series. The implementation of a generic source of identical design in all TPSs using MBDCAs is an important step toward supporting univocal commissioning procedures and direct comparisons between TPSs.

BrachyGuide: a brachytherapy-dedicated DICOM RT viewer and interface to Monte Carlo simulation software.

This work presents BrachyGuide, a brachytherapy-dedicated software tool for the automatic preparation of input files for Monte Carlo simulation from treatment plans exported in DICOM RT format, and results of calculations performed for its benchmarking. Three plans were prepared using two computational models, the image series of a water sphere and a multicatheter breast brachytherapy patient, for each of two commercially available treatment planning systems: BrachyVision and Oncentra Brachy. One plan involved a single source dwell position of an 192Ir HDR source (VS2000 or mHDR-v2) at the center of the water sphere using the TG43 algorithm, and the other two corresponded to the TG43 and advanced dose calculation algorithm for the multicatheter breast brachytherapy patient. Monte Carlo input files were prepared using BrachyGuide and simulations were performed with MCNP v.6.1. For the TG43 patient plans, the Monte Carlo computational model was manually edited in the prepared input files to resemble TG43 dosimetry assumptions. Hence all DICOM RT dose exports were equivalent to corresponding simulation results and their comparison was used for benchmarking the use of BrachyGuide. Monte Carlo simulation results and corresponding DICOM RT dose exports agree within type A uncertainties in the majority of points in the computational models. Treatment planning system, algorithm, and source specific differences greater than type A uncertainties were also observed, but these were explained by treatment planning system-related issues and other sources of type B uncertainty. These differences have to be taken into account in commissioning procedures of brachytherapy dosimetry algorithms. BrachyGuide is accurate and effective for use in the preparation of commissioning tests for new brachytherapy dosimetry algorithms as a user-oriented commissioning tool and the expedition of retrospective patient cohort studies of dosimetry planning.

Experimental determination of the Task Group-43 dosimetric parameters of the new I25.S17plus (125)I brachytherapy source.

PURPOSE: To present experimental dosimetry results for the new IsoSeed I25.S17plus (125)I brachytherapy source, in fulfillment of the American Association of Physicists in Medicine recommendation for, at least one, experimental dosimetry characterization of new low-energy seeds before their clinical implementation. METHODS AND MATERIALS: A batch of 100 LiF thermoluminescent dosimeter (TLD)-100 microcubes was used for the experimental determination of the dose-rate constant, radial dose, and anisotropy functions, in irradiations performed using two Solid Water phantoms. Monte Carlo (MC) simulations were used to determine appropriate correction factors that account for the use of Solid Water as a phantom material instead of liquid water and for the different energy response of the TLD dosimeters in the experimental (125)I photon energies relative to the 6 MV x-ray photon beam used for the TLD calibration. Measurements were performed for four I25.S17plus seeds; one with direct traceability of air-kerma strength calibration to National Institute of Standards and Technology and three with secondary National Institute of Standards and Technology traceability. RESULTS: A mean dose-rate constant, Λ, of 0.956 ± 0.043 cGy h(-1) U(-1) was experimentally determined for the I25.S17plus source, which agrees within uncertainties with the MC result of 0.925 ± 0.013 cGy h(-1) U(-1) calculated independently for the same seed model in a previous study. Agreement was also observed between the measured and the MC-calculated radial dose and anisotropy function values. CONCLUSIONS: Experimental dosimetry results for the I25.S17plus (125)I source verify corresponding independent MC results in the form of Task Group-43 dosimetry parameters. The latter are found in agreement within uncertainties with sources of similar design incorporating a silver marker, such as the Oncura OncoSeed Model 6711.

New (125)I brachytherapy source IsoSeed I25.S17plus: Monte Carlo dosimetry simulation and comparison to sources of similar design.

PURPOSE: To determine the relative dose rate distribution around the new (125)I brachytherapy source IsoSeed I25.S17plus and report results in a form suitable for clinical use. Results for the new source are also compared to corresponding results for other commercially available (125)I sources of similar design. MATERIAL AND METHODS: Monte Carlo simulations were performed using the MCNP5 v.1.6 general purpose code. The model of the new source was prepared from information provided by the manufacturer and verified by imaging a sample of ten non-radioactive sources. Corresponding simulations were also performed for the 6711 (125)I brachytherapy source, using updated geometric information presented recently in the literature. The uncertainty of the dose distribution around the new source, as well as the dosimetric quantities derived from it according to the Task Group 43 formalism, were determined from the standard error of the mean of simulations for a sample of fifty source models. These source models were prepared by randomly selecting values of geometric parameters from uniform distributions defined by manufacturer stated tolerances. RESULTS AND CONCLUSIONS: Results are presented in the form of the quantities defined in the update of the Task Group 43 report, as well as a relative dose rate table in Cartesian coordinates. The dose rate distribution of the new source is comparable to that of sources of similar design (IsoSeed I25.S17, Oncoseed 6711, SelectSeed 130.002, Advantage IAI-125A, I-Seed AgX100, Thinseed 9011). Noticeable differences were observed only for the IsoSeed I25.S06 and Best 2301 sources.

The effect of finite patient dimensions and tissue inhomogeneities on dosimetry planning of 192Ir HDR breast brachytherapy: a Monte Carlo dose verification study.

PURPOSE: To evaluate the accuracy of clinical dosimetry planning using commercially available treatment planning systems in (192)Ir high-dose-rate (HDR) breast brachytherapy, with emphasis on skin dose, in view of potential uncertainties owing to the patient finite dimensions and the presence of the lung. METHODS AND MATERIALS: A patient-equivalent mathematical phantom was constructed on the basis of the patient computed tomography scan used in the clinical treatment planning procedure. The actual treatment plan delivered to the patient, involving an implant of six plastic catheters and 26 programmed source dwell positions, was simulated by means of the Monte Carlo method. Results are compared with corresponding dose calculations of a commercially available treatment planning system in the form of prescribed dose percentage isodose contours and cumulative dose-volume histograms. RESULTS: The comparison of Monte Carlo results and treatment planning system calculations revealed that all percentage isodose contours greater than 60% of the prescribed dose are not affected by the finite breast dimensions or the presence of the lung. Treatment planning system calculations overestimate dose in the lung as well as lower isodose contours at points lying both close to the breast or lung surface and relatively away from the implant. In particular, skin dose is overestimated by 5% in the central breast region and within 10% at all other points. CONCLUSIONS: Dose-volume histogram and all other relevant planning quality indices for the planning target volume calculated by the treatment planning system are credible. Skin and lung dose calculations by the treatment planning system can be thought of as a conservative approach in view of the reported dose overestimation.

On the dosimetric accuracy of a Sievert integration model in the proximity of 192Ir HDR sources.

PURPOSE: To investigate the efficacy of a Sievert integration model in dosimetry close to 192Ir high-dose-rate brachytherapy sources and validate its accuracy and potential to resolve dosimetric differences between these sources in the cm and mm distance ranges relevant to interstitial and intravascular brachytherapy applications, respectively. METHODS AND MATERIALS: The dosimetric quantities of the generalized Task Group 43 formalism, as well as dose rate profiles in polar and Cartesian coordinates, are calculated, and results are compared to corresponding Monte Carlo data in the literature. RESULTS: Sievert calculations were found in excellent agreement with corresponding Monte Carlo published results. Dose rate polar angle profiles in the cm distance range depended significantly on corresponding anisotropy function data, whereas in the mm distance range, dose rate polar angle profiles are governed by the corresponding geometry function profiles, because anisotropy proved insignificant. Radial dose functions of the sources were found comparable. A simple equation for the calculation of the dose rate constant of the sources within clinically acceptable accuracy is provided. CONCLUSIONS: The particular Sievert model proved capable of resolving dosimetric differences of the sources and provides results within clinical accuracy. Therefore, it constitutes a useful tool for dosimetry in clinical practice and especially in intravascular applications, where there is currently a lack of available dosimetric data.